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Introduction 


b-lsonucleotide (isoNA) incorporation around 
cleavage site of passenger strand promotes the 
vibration of Ago2-PAZ domain and enhances 
in vitro potency of siRNA 


Ye Huang,t®* Miao Tian, t° Yichao Zhang,” Gang Sheng, Zhuo Chen,? Yuan Ma,? 
Yue Chen,? Yihong Peng,“ Yi-Lei Zhao,*? Yanli Wang,‘ Lihe Zhang? and Zhenjun Yang** 


It has been demonstrated that passenger strand cleavage is important for the activation of RNA-induced 
silencing complex (RISC), which is a crucial step for siRNA-mediated gene silencing. Herein, we report 
that isonucleotide (isoNA) modification around the cleavage site of the passenger strand would affect the 
in vitro potency of modified siRNAs by altering the motion pattern of the Ago2-PAZ domain. According to 
western blotting, q-PCR and antiviral test results, we proved that p-isonucleotide (isoNA) modification at 
the position 8 of the passenger strand (siMek1-SO8D), which is adjacent to the cleavage site, markedly 
improved the in vitro potency of the modified siRNA, whereas siRNAs with D-isoNA incorporation at posi- 
tion 9 (siMek1-SO9D) or L-isoNA incorporation at positions 8 and 9 (siMek1-SO8L, siMek1-SO9L) displayed 
lower activity compared to native siRNA. Kinetics evaluation of passenger strand cleavage induced by 
T. thermophilus Ago (Tt-Ago) showed that p-isoNA modification at position 8 of the passenger strand had 
no significant influence on the cleavage rate, but L-isoNA modification at position 8 slowed the cleavage 
rate markedly. Moreover, the results of molecular dynamics simulations showed that D-isoNA modification 
at position 8 affected the open-close motion of the PAZ domain in the Ago/siRNA complex, which may 
promote the loading of RISC and release of a passenger strand cleavage product, and consequently 
accelerate the activation of RISC and enhance silencing activity. However, p-isoNA modification at 
position 9 or L-isoNA modification at position 8 or 9 exerted opposite influences on the motion of the 
Ago-PAZ domain. 


(RISC) is the pivotal step in the RNAi pathway. In humans, a 
siRNA duplex with both guide strand and passenger strand is 


Since the discovery of the RNA interference (RNAi) effect in 
Caenorhabditis elegans,‘ double-stranded small interfering 
RNA (siRNA) molecules have been proved to be able to exert a 
potent gene knockdown effect in mammalian cells” and other 
eukaryotes?” Assembly of the RNA-induced silencing complex 


“State Key Laboratory of Natural and Biomimetic Drugs, School of Pharmaceutical 
Sciences, Peking University, Beijing 100191, China. E-mail: yangzj@bjmu.edu.cn 
State Key Laboratory of Microbial Metabolism, School of Life Sciences and 
Biotechnology, Shanghai Jiao Tong University, Shanghai, China. 

E-mail: yileizhao@sjtu.edu.cn 

“Laboratory of Non-Coding RNA, Institute of Biophysics, Chinese Academy of 
Sciences, Beijing, China 

“Department of Microbiology, School of Basic Medical Sciences, Peking University, 
Beijing, China 

“Department of New Drug Research and Development, Institute of Materia Medica, 
Peking Union Medical College & Chinese Academy of Medical Sciences, Beijing, 
China 

+ These authors contributed equally to this paper. 


This journal is © The Royal Society of Chemistry 2015 


first recognized by the RLC complex, which contains Dicer and 
the RNA-binding protein TRBP, to form an immature RISC 
complex.’ Then, the passenger strand is nicked between posi- 
tions 9 and 10, counted from the 5’-terminal, by the Argonaute 
2 (Ago2) protein in the same way as target mRNA is cleaved,*” 
followed by being degraded by C3PO to finally activate the 
RISC.'° Active RISC recognizes and cleaves complementary 
target RNA to execute the silencing process. 

A dramatic reduction in the cleavage efficiency of target 
RNA was found when chemical modification occurred around 
the nicked site of the passenger strand. These chemical modi- 
fications may limit the formation of RISC and then impair the 
silencing activity of siRNA by blocking the cleavage of the pas- 
senger strand.*"' For example, replacing the phosphodiester 
linkage between position 9 and position 10 with phos- 
phorothioate, modifying position 9 with 2’-O-methylribose, 
and introducing 4-nt mismatches from position 9 to position 
11, all caused a great reduction in the cleavage efficiency of 


Org. Biomol. Chem., 2015, 13, 10825-10833 | 10825 


Paper 


complementary target RNA. However, the same modifications 
that were positioned 1- or 2-nt downstream of the cleavage site 
or away from the central position did not affect the cleavage of 
target RNA."* On the other hand, when thermally destabilizing 
modifications were placed around the cleavage site of the pas- 
senger strand, an obvious enhancement in silencing activity 
was found.” This unusual phenomenon may be due to the 
fact that thermally destabilizing modifications change the 
regional conformation around the active site and assist 
the cleavage-related residues of Ago2 protein to interact with 
the phosphodiester bond between position 9 and position 
10 of the passenger strand. All these reports indicate that 
alteration of the conformation around the cleavage site of the 
passenger strand may affect the silencing activity. 
Isonucleotide (isoNA) (Fig. 1A) is a novel type of nucleotide 
analogue in which the nucleobase is moved to another posi- 
tion of ribose other than C-1'."** In our previous study, we 
found that isonucleotide modification could alter the regional 
conformation around the incorporation site when modified 
DNA or RNA oligonucleotides were hybridized with a com- 
plementary DNA/RNA strand. For example, t-isoNA modification 
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at the 3’-terminus of DNA oligonucleotides created a con- 
formation that blocked exonuclease recognition and then 
increased their exonuclease-resistant ability, whereas L-isoNA 
modification in the middle of oligonucleotides altered the 
regional conformation into a form that is more suitable for 
RNase H recognition and then promoted the degradation of 
complementary RNA.*°*” Then, we reported that p-isoNA and 
L-isoNA incorporation could introduce regional conformation 
alteration into a siRNA duplex and affect the physicochemical 
and biological properties at varying degrees. Therefore, 
p-/L-isoNA could be applied as a pair of chemical tools to inves- 
tigate the systematic relationship between conformational 
alteration around every single nucleotide of siRNA and their 
biological potency. 

Through thermal denaturation and dynamic simulation 
studies, we have demonstrated that p-isoNA and L-isoNA modi- 
fication at positions 8 and 9 of the passenger strand would 
affect the original compact duplex by introducing flexible con- 
formations of different degrees. Therefore, in this study, 
further silencing activity, antivirus activity and dynamic simu- 
lation studies were applied to thoroughly explore the influence 
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Fig. 1 Silencing activity affected by D-/L-isoNA modification at positions 8 and 9 of passenger strand. (A) Structures of p-isonucleotide (D-isoNA, 
left) and L-isonucleotide (L-isoNA, right). (B) Sequences of siMek1s modified with p-isoNA or L-isoNA in the passenger strand. The blue dotted lines 
represent the predicted cleavage sites. (C) Quantitative PCR results of mek1 mRNA. HEK293 cells were co-transfected with siMek1 and siMekis 
modified with D-isoNA or L-isoNA in the passenger strand at 30 nM of each. Mek1 mRNA was measured 24 h later by quantitative PCR. Mek1 mRNA 
levels were normalized to B-actin mRNA, which was used as a loading control. The uninhibited normalized mek1 mRNA level was set to 100%. The 
data are the average of three independent experiments and error bars denote the standard deviations. (D) Western blotting analysis of MEK1 protein. 
HEK293 cells were co-transfected with siMek1 and siMek1s modified with p-isoNA or L-isoNA in the passenger strand at 30 nM of each. Protein 
lysates were prepared for 36 h and blotted with MEK1 antibodies. Quantification of mek1 bands was normalized by densitometric scanning of actin 
using the Band Leader software (version 3.0). The relative intensity of MEK1 on the y-axis is represented by the percentage of the intensity of mock 
transfection (100%). The data are the average of three independent experiments and error bars denote the standard deviations. (E) Western blotting 
analysis of MEK1 protein when HEK293 cells were co-transfected with siRNAs at 5 nM, 10 nM, and 30 nM. 
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of regional conformation alteration caused by p-/L-isoNA modi- 
fication on the process of passenger strand cleavage and the 
silencing potency of modified siRNAs. The siRNA sequence 
that targets the MEK1 kinase mRNA, which is crucial for the 
replication process of various common viruses, was chosen 
for the further evaluation of silencing activity and antivirus 
activity in vitro (Fig. 1B).'?° 


Results 


Silencing activity affected by p-/L-isoNA incorporation at 
positions 8 and 9 of passenger strand 


Quantitative PCR and western blotting were utilized to 
investigate the silencing activity of the synthesized siRNAs with 
p-/L-isoNA modification at passenger strand at RNA level and 
protein level, respectively, in a HEK293 cell line. According to 
the results of qPCR evaluation, when position 8 on the passen- 
ger strand was replaced with p-isoNA (siMek1-SO8D), greater 
silencing potency than the native sequence was found; on 
the contrary, t-isoNA at the same position 8 (siMek1-S08L) 
impaired the silencing potency. However, at the position 9, 
both p-isoNA and t-isoNA exhibited decreased silencing 
activity, although the degree of decreasement was not as 
obvious as siMek1-SO8L (Fig. 1C). This same phenomenon 
was also found at protein level, namely, SO8D exhibited the 
strongest knockdown effect, whereas SO8L exhibited worse 
effect than the native sequence. p-/t-isoU modification at the 
position 9 showed decreased silencing activity compared to 
the native sequences and the difference was more obvious at 
protein level (Fig. 1D). The difference between the silencing 
activities of siMek1-SO8D, siMek1-SO8L and native siMek1 were 
nearly the same at three different concentrations of siRNA, 
namely, 5 nM, 10 nM, and 30 nM. This indicated that siMek1- 
SO8D provided the best silencing activity regardless of the con- 
centration, although the degree of enhancement was not 
obvious when the concentration of siRNA was too low, such as 
5 nM (Fig. 1E). In a word, p-isoNA modification at the position 
8 of the passenger strand was able to enhance the in vitro 
potency of siRNA. 


Influence of modification with isoNA on the kinetics of 
passenger strand cleavage induced by T. thermophilus Ago 
(Tt-Ago) 


According to the results of the evaluation of silencing activity, 
modification with p-isoNA at the position 8 of the passenger 
strand may improve the silencing effect of siRNA. Conversely, 
other modification strategies around this area, such as the 
incorporation of L-isoNA at the same position or p-/L-isoNA at 
the adjacent position 9, would impair the activity. In previous 
reports, the chemical modifications that disturb cleavage of 
the passenger strand were found to be able to impair the silen- 
cing activity of siRNA by limiting the formation of RISC.* This 
reminds us that modification with isoNA around the cleavage 
site of the passenger strand may also modulate silencing 
activity by directly affecting the efficiency of the cleavage 
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process or other steps. Therefore, the influence of modifi- 
cation with isoNA on the kinetics of passenger strand cleavage 
induced by T. thermophilus Ago (Tt-Ago) was evaluated. Pre- 
vious reports studied the cleavage kinetics of the passenger 
strand with a 5’-°’P-radiolabeled RNA oligonucleotide.’ They 
detected the cleavage fragment of the passenger strand via 
autoradiography imaging of *’P. However, the labelling process 
using °P is complicated, time-consuming and dangerous. 
Herein, we detected the cleavage fragment with 3’-terminal-Cy3- 
labelled passenger RNA instead of a 5’-*’P-radiolabeled RNA 
oligonucleotide. To avoid the influence caused by the steric 
hindrance of the Cy3 group on RISC activation, four extra 
uridine nucleotides were added at the 3'-terminus (Fig. 2A). 

According to the results of cleavage kinetics (Fig. 2B), when 
p-isoNA and L-isoNA were incorporated at positions 8 and 9 of 
the passenger strand, neither of them disturbed the cleavage 
process dramatically. The cleavage of all five 3’-Cy3-labelled 
RNA strands, Cy3-MKS-23, Cy3-MKS-2308D, Cy3-MKS-2308L, 
Cy3-MKS-2309D, and Cy3-MKS-2309L, started after 30 min and 
finished after 120 min. However, according to the fluorescence 
intensity at 60 min, we found that modification with p-isoNA 
affected the cleavage process lesser than modification with 
L-isoNA. Cy3-MKS-2308D and Cy3-MKS-2309D had nearly the 
same cleavage rate as Cy3-MKS-23, but the cleavage of 
Cy3-MKS-2308L and Cy3-MKS-2309L was slower. This indi- 
cated that the low silencing activity of SO8L and SO9L may be 
due to the disturbance caused by t-isoNA incorporation to the 
interaction between the PIWI domain of Ago and the cleavage 
site of passenger RNA. In addition, the stronger silencing 
activity of SO8D could not be simply explained through the 
cleavage kinetics experiment. 


Stability of the Ago/siRNA complexes and change in motion of 
Ago functional domains 


Then, the mechanism of this difference was investigated by 
molecular dynamics simulations. RMSF of five simulated 
models of a complex of Ago2 and siRNAs are shown in Fig. 3A. 
In the Ago_Native model with no modifications, there are two 
loops that are obviously flexible. As shown in Fig. 3B, the two 
loops belong to two different functional domains, PAZ and 
PIWI, and the PAZ loop is more flexible than the PIWI loop. 
After isoNA modifications, the flexibility of the two loops 
decreased, especially for the PAZ loop. Ago_SO08D model 
exhibited the smallest decrease, whereas Ago_SO8L displayed 
the largest decrease. For the siRNA part, there was no detect- 
able change in flexibility. 

The four different domains of the Ago protein have crucial 
functions in the cleavage process. We found that the flexibility 
of the PAZ and PIWI domains changed after isoNA modifi- 
cation, but there is no information regarding how the motions 
of the domains vary. Principle component analysis (PCA) was 
performed on the Ago/siRNA models to reveal the significant 
transformations of the functional domains. From the PC1 in 
Ago_Native model (Fig. 3D), we found an obvious open-close 
motion for the PAZ loop, which is in accordance with the 
RMSF result and is also found in other molecular dynamics 
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A Predicted Cleavage Site 
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Influence of isoNA modification on the kinetics of passenger strand cleavage induced by T. thermophilus Ago (Tt-Ago). (A) Sequences of 3’- 


Cy3-labeled p-/L-isoNA-modified siRNA. (B) Kinetics of passenger strand cleavage induced by T. thermophilus Ago (Tt-Ago). Reaction buffer: 10 mM 
Hepes (pH = 7.5), 100 mM NaCl, 5 mM MnCl; 5 pM Tt-Ago and 5 uM guide DNA strand, incubation at 42 °C for 30 min; then, 3 M passenger RNA 
strand was added and the incubation temperature was increased to 70 °C for cleavage reaction; the cleavage mixtures were sampled at 0 min, 
15 min, 30 min, 60 min, 120 min, and 180 min, and separated on a 20% PAGE gel with cleavage fragments and visualized under fluorescent light. 
The black arrow points at the intact passenger RNA substrate and the red arrow points at the 3'-cleavage passenger RNA fragment. 


simulations. The open-close motion of the PAZ loop is thought 
to be important for the propagation of the guide and target 
strand and tends to be significantly different, depending on 
whether the target strand participates or not. It is conceivable 
that the open-close motion of the PAZ loop is essential for the 
uptake of siRNA duplex, and the propagation of the guide and 
target strand may also determine the release of passenger 
strand cleavage products. For the Ago_SO08D model, similar 
motions were found in PC2 and PC3, which means that the 
PAZ loop moves with smaller eigenvalues and thereby higher 
frequencies. To confirm the open-close motion, we determined 
the distance between the most flexible residue in the PAZ loop 
and the active site. The open-close ranges vary from 30 A to 
57 A along this indicator (Fig. 3C) for the Ago_Native model. A 
similar motion was also observed for the Ago_SO8D model but 
with a slightly smaller cycle. The other three models do not 
exhibit similar motions and the distances change little. On the 
other hand, there is no obvious motion for the PIWI loop 
around 490, which has high RMSF; this means that it behaves 
randomly. 


Change in binding affinity after isoNA modifications 


The binding free energy of RISC/ssRNA indicates the release 
obstacle of the product after cleavage of the passenger strand. 
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After isoNA modifications, the binding free energy decreases 
for Ago_SO8D, whereas those of the other three modifications 
(Ago_SO8L, Ago_SO9D and Ago_SO9L) increase (Table 1). The 
activation of RISC may become easier for Ago_SO8D due to the 
easier discarding of the passenger strand, whereas the other 
three isoNA modification strategies prevent the formation of 
active RISC. The difference may be related to the motion of the 
PAZ loop, where the channel with SO8D involved closed 
complex more than the native complex (Fig. 3) and then 
inhibited the interaction between S08D and Ago2 protein. 


Antiviral effects of p-/L-isoNA incorporation at position 8 of 
passenger strand 


Then, the difference between the silencing activities of 8D and 
8L was further investigated via the downstream antiviral effect. 
Human herpes simplex virus type 2 (HSV2) belongs to the 
alpha herpesvirus subfamily and infection with HSV2 pro- 
duces high morbidity and exhibits broad prevalence, long 
latency, and difficult curability. To date, no preventive vaccines 
or effective therapeutic drugs for recurrent infections are avail- 
able.** Enterovirus 71 (EV71) is a positive single-stranded RNA 
virus belonging to the family of Picornaviridae, which is 
known to cause hand, foot and mouth disease (HFMD) in 
young children and can lead to the development of severe 
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Fig. 3 Stability of the Ago/siRNA complexes. (A) RMSF of each residue in the Ago/siRNA complexes for the last 20 ns; (B) Ago/siRNA complex of 
Ago_Native model. The red loop refers to the first RMSF peak that appears around residue 220, and the blue loop refers to the second RMSF peak 
that appears around residue 490; (C) distances between the most flexible residue and the active site of each model during the simulation; (D) the 
first three principal components (PCs) of Ago_Native and Ago_SO8D models during the simulation period. 


Table 1 The binding free energies of Ago/siRNA and RISC/ssRNA for 
each of the models (kcal mol?) 


Ago/siRNA RISC/ssRNA 
Ago_Native —587.7944 + 25.5832 —254.2087 + 20.9272 
Ago_S08D —567.7166 + 24.0868 —246.0320 + 17.8012 
Ago_SO8L —603.6906 + 32.4930 —282.1824 + 24.5516 
Ago_S09D —610.9651 + 22.4295 —277.2535 + 15.9013 
Ago_SO9L —598.7532 + 25.6560 —278.9986 + 19.8106 


neurological diseases in some patients, especially children. To 
date, the treatment and control of EV71 infection are only 
symptomatic due to the lack of effective medications and un- 
availability of a prophylactic vaccine.” The inhibition of MEK 
has been studied as the most effective way of investigating the 
importance of the MAPK cascade in various cellular physio- 
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logical and pathological processes, including HSV2 and EV71 
infection. Herein, we studied the anti-HSV2 and anti-EV71 
effects of siMek1 with modification by p/t-isoA at position 8 of 
the passenger strand by western blotting, virus titration and 
morphological analysis. To further determine the anti-viral 
effect of modified siMek1s on the activation of ERK by HSV2 
or EV71, HEK293 cells or RD cells were transfected with 
siRNAs and then infected with HSV2 or EV71, respectively. As 
shown in Fig. 4A and B, siMekis effectively silenced the 
expression of MEK1 and viral gB protein. And, SO8D had a 
better antiviral effect than either the native sequence or SO8L, 
which is in accordance with the former results of silencing 
activity evaluation. 

As one of the most important indices of viral replication, 
viral titer reflects the release of the live virus in the super- 
natant of cells infected with the virus. To investigate the influ- 
ence of p-/t-isoNA-modified siRNAs in HSV2 replication, 
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Fig. 4 Antiviral effects of p-/L-isoNA incorporation at position 8 of passenger strand. (A) and (B) Anti- HSV2 and anti-EV71 results of p-/L-modified 
siMekis. The experiments were performed as described in Fig. 1 except for virus infection with an MOI of 1.5. Cell lysates that were collected at 24 h 
p.i. were assayed for viral protein expression by western blot with monoclonal antibodies specific to HSV2 gB protein or EV71 VP1 protein. Quantifi- 
cation of the gB band and VP1 band was normalized by densitometric scanning of actin using the Band Leader software (version 3.0). The relative 
intensity of gB and VP1 on the y-axis is represented by the percentage of the intensity of mock transfection (100%). The data are the average of three 
independent experiments and error bars denote the standard deviations. (C) Viral titer of HSV2 and EV71 after incubation with p-/L-isoNA-modified 
siMekis. HEK293 cells were infected with HSV2 at an MOI of 1.5. Cells and medium were collected separately at different time points. Virus titers in 
the medium were measured by plaque assay. Data shown are means + standard deviations (n = 3). The titers on the y-axis represent logarithmic 
values. (D) CPE results with HSV2 and EV71 infection after incubation with p-/L-isoNA-modified siMekis. Cells were transfected with modified 
siRNAs (30 nmol of each) and infected at 44 h post-transfection with 1.5 MOI of HSV2. The effects of modified siRNAs on CPE were examined from 


8 hp.i. and images were obtained at 24 h p.i. 


HEK293 cells transfected with siMekis were infected with 
HSV2, whereas RD cells were infected with EV71 and com- 
pared with control cells that were infected with HSV2 or EV71 
alone, and then the virus production was examined. As shown 
in Fig. 4C, in the cells that were transfected with 30 nM 
siRNAs, HSV2 titers were markedly inhibited as compared to 
control cells that were infected with HSV2 or EV71 alone. 
Although not as obviously as HSV2, EV71 titers were also 
inhibited by siMekis. Moreover, an enhanced anti-viral effect 
can be seen in SO8D compared with unmodified siRNA and 
SO8L. 

To determine the effects of siMek1s on the cytopathic effect 
(CPE) of HSV2 and EV71, HEK293 cells or RD cells that were 
transfected with siMekis were infected with HSV2 and EV71 
and then monitored by phase-contrast microscopy every 12 h 
for 3 days p.i. In this experiment, mock infection and viral 
infection were used as controls. As shown in Fig. 4D, at 24 h, 
p.i. cells that were transfected with 30 nM siMek1s manifested 
considerably less CPE as compared to virus and siMut1 con- 
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trols. Notably, SO8D also exhibited less CPE than the others. 
The virus titer that was observed in the supernatant further 
confirmed that modification with p-isoNA displayed some 
potency-enhancing effect at passenger strand position 8, 
which could be applied in anti-HSV2 and EV71 therapy. 


Discussion 


Generally, during the process of RNA interference, the cleavage 
of the passenger strand caused by Ago2 protein occurs 
between nucleotide 9 and nucleotide 10. Chemical modifi- 
cations and mismatches located at these two sites are able to 
interfere with the cleavage process by disturbing the approach 
of and interaction between cleavage-related residues in the 
Ago2-PAZ domain and the phosphodiester bond between 
nucleotide 9 and nucleotide 10. As a result, the silencing 
activity of modified siRNA would also be reduced."* In this 
report, we also found that irrespective of the incorporation of 
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p-isoNA or t-isoNA at position 9, the silencing efficiency of 
modified siRNA was reduced markedly, which is in consistency 
with other modification strategies. However, chemical modifi- 
cations at other positions usually have no obvious influence 
on the cleavage process of the passenger strand.** However, 
we found that p-isoNA and t-isoNA modification next to the 
cleavage-related site could also affect the cleavage process. 
Interestingly, p-isoNA and t-isoNA modification introduce con- 
siderably opposite influences on the silencing efficiency when 
the modification occurs at position 8 of the passenger strand. 

According to the cleavage kinetics experiment, the passen- 
ger RNA strand with modification by t-isoNA was cleaved con- 
siderably slower than the unmodified passenger strand at the 
same position. Based on our previous reports,’® L-isoNA modi- 
fication may introduce dramatic alterations in regional confor- 
mation not only at the incorporation site, but also around the 
incorporation site. Therefore, when t-isoNA modification was 
placed next to the nucleotide 9, it also disturbed the tight 
interaction between the cleavage-related residues and the 
corresponding phosphodiester bond. Although in the direct 
cleavage kinetics test, the passenger RNA strand with p-isoNA 
modification at position 8 exhibited no difference from the 
unmodified passenger strand; they exerted quite different 
influences on the vibration frequency of the Ago2-PAZ domain 
in the molecular dynamics simulation. This may also be due 
to the slight alteration in regional conformation around the 
insertion site caused by modification with p-isoNA. However, 
study on the precise mechanism of the enhanced vibration fre- 
quency of the Ago2-PAZ domain is still ongoing, and to the 
best of our knowledge, no modification strategy has enhanced 
the silencing activity by increasing the movement frequency of 
the PAZ domain. Although the degree of activity enhancement 
was not quite dramatic, it provided us with more accurate 
knowledge about the RNA interference process. Our study also 
reminds us that p-isoNA and t-isoNA, as a pair of chemical 
tools, could exert quite opposite influences on the interaction 
between proteins and nucleotides, and they can be applied 
widely in mechanism studies in the future. 


Experimental 
Solid-phase synthesis of RNA oligonucleotides 


Synthesis of RNA oligonucleotides was carried out at the 
1 pmol scale using an Applied Biosystems 394 DNA synthesizer 
according to regular phosphoramidite chemistry. 2'-OTBDMS- 
protected pD- and t-isoNA phosphoramidite monomers were 
synthesized as previously reported? and then the p-/L-isoNA 
phosphoramidite monomers were inserted singly into the 
siRNA sequence at the modified position targeting MEK1 
(Fig. 1B). Furthermore, an extended coupling time of 900 
seconds was used for p-/L-isoNA phosphoramidites instead of 
the standard coupling time of 600 seconds that is used for the 
four standard 2’-OTBDMS-protected RNA phosphoramidites 
(A®”, C^, G^ and U) due to the steric effect of isoNA phosphor- 
amidites. Cleavage and deprotection of the oligomers were 
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performed in concentrated ammonia at 55 °C for 24 h and the 
TBDMS group was removed by treatment with Et;N/HF in 
anhydrous DMSO. The crude product was purified by anion- 
exchange high-performance liquid chromatography (Dionex 
DNAPac PA200, 9 x 250 mm) using a linear gradient of 12.5%- 
40% eluent A for 40 min. A solution of 0.02 M Tris-HCI1O, in 
10% CH3CN, pH = 8, was used as eluent B, and 0.4 M NaClO, 
in eluent B was used as eluent A. Then, the purified oligonu- 
cleotides were desalted by a Sephadex G25 column. The purity 
and identity of the oligonucleotides were confirmed by ion- 
exchange chromatography and MALDI-TOF MS, respectively, 
and all single strands were purified by HPLC to >90% purity. 
Finally, the pure oligonucleotides were lyophilized and stored 
at —40 °C. 


Antibodies 


Antibodies were purchased from Cell Signaling Technology 
(anti-MEK1), EastCoast Bio (anti-HSV gB monoclonal 
antibody), Abcam (polyclonal anti-EV71 VP1) and Santa Cruz 
Biotechnology (anti-B-actin). 


Cell culture 


Rhabdomyosarcoma (RD) cells and human embryonic kidney 
(HEK) 293 cells were grown in Dulbecco’s modified Eagle’s 
medium (DMEM, GIBCO) supplemented with 2% or 10% heat- 
inactivated fetal bovine serum (FBS, GIBCO) at 37 °C in a 5% 
CO, atmosphere in a humidified incubator. 


Virus infection 


Monolayers of RD or HEK293 cells at 70%-80% confluency 
were starved and infected with strain 333 of HSV2 and wild- 
type EV71 (EV71-BC08 strain) at various multiplicities of infec- 
tion (MOI). After 1 h at 37 °C, the inoculum was removed and 
incubation was continued in DMEM containing 2% FBS at 
37 °C for viral growth. Supernatants that were collected at 
different time points post-infection (p.i.) were clarified by cen- 
trifugation (2000g). Virus titers were determined by plaque 
assays carried out in triplicate. Viruses were inactivated by 
UV-irradiation for 30 min using a Phillips TUV 30 W bulb at a 
distance of 17 cm with occasional agitation, and then blind- 
passaged for three generations to confirm the loss of replica- 
tive capacity. 


siRNA and transient transfection 


For siRNA transfection, the INTERFER in reagent (Polyplus- 
transfection) was used according to the manufacturer’s 
instructions. RD cells or HEK293 cells were grown in 12-well 
plates to 40% confluency. Cells were then transfected with 
30 nM specific siRNAs and incubated for the indicated 
number of hours. The medium was then replaced by DMEM 
containing 1% FBS and the cells were starved for 8 h. 


Morphological analysis 


HEK293 cells were examined following infection by HSV2 and 
EV71 at 12 h intervals post-infection for cytopathic effects 
(CPE) using phase-contrast microscopy. 
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Western blotting analysis 


Cells were lysed in a buffer containing 20 mM Hepes (pH 7.4), 
100 mM NaCl, 5 mM EDTA (pH 7.4), 1 mM Na3VO,, 30 mM 
NaF, 5% glycerol, 0.1% SDS, 1% Triton X-100, 10 mM p-nitro- 
phenyl phosphate, and 1 mM glycerophosphate, supplemented 
with complete protease inhibitors (Roche). Cell lysates were 
obtained by centrifugation at 4 °C and 13 000 rpm and protein 
concentrations were determined using the Bicinchoninic Acid 
Protein Assay Kit (Pierce). Proteins were resolved onto sodium 
dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) 
gels and transferred to PVDF membranes (Millipore). The 
membrane was blotted with specific primary antibodies, as 
indicated in the figure legends, followed by incubation with 
secondary antibodies conjugated with horseradish peroxidase 
and developed with an enhanced chemiluminescence sub- 
strate (ECL) or SuperSignal West Femto Maximum Sensitivity 
Substrate (Pierce). 


Real-time quantitative PCR (qPCR) 


Viral RNAs from supernatants were prepared using Trizol 
reagent (Invitrogen). Then, 1 pg total RNA was reverse tran- 
scribed into cDNA using the RevertAid First Strand cDNA syn- 
thesis kit (Thermo Scientific). A quantitation standard curve 
was obtained using seven 10-fold serial dilutions of recombi- 
nant plasmid standard DNA, which contained sequences 
derived from the gB gene of HSV2 and the VP1 gene of EV71, 
varying from 1 x 10° to 1 x 10° copies. Real-time quantitative 
PCR (qPCR) was performed by targeting the gB gene and 
VP1 gene in 96-well plates using the Roche LightCycler 480 
system. Each 20 pL reaction mixture contained 10 pL 2x SYBR 
Green Supermix, 20 nM target gene primer mix, and 20-50 ng 
cDNA template. The PCR reaction was set up as follows: initial 
denaturation step at 95 °C for 10 min, followed by 40 cycles of 
30 s at 94 °C, at 55 °C for 30 s and at 72 °C for 30 s. Quantified 
results for the experimental samples were extrapolated from 
the standard curve, with all experimental samples being run in 
triplicate. 


Molecular structures of Ago/siRNA complex 


Eukaryotic Ago2 protein crystal structures were as reported,*” 
whereas siRNA double helixes are not involved in the struc- 
tures. Owing to the high structural similarity of the Ago family, 
Ago protein from T. thermophilus was used as a substitute. 
There are many T. thermophilus Ago proteins in the Protein 
Data Bank such as 3HJF, 3HK2, 3HM9, 3HVR, 3HO1” and 
3F73.”* Ago protein from 3HVR was used because there are no 
mutations and mismatches and there are two Mg” ions in the 
active site. The total nucleotide part from 3HK2 was used after 
the overlap of the Ago proteins in 3HVR and 3HK2, in which 
the RMSD of the same guide and target nucleotides (27 nt) 
between 3HVR and 3HK2 is 1.54 A. Finally, some processes 
were carried out to make the siRNA available for the following 
molecular simulation, which included changing the DNA 
guide strand to an RNA guide strand, mutating the nucleotide 
sequence and modifying a specific nucleotide to an isonucleo- 
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tide according to Fig. 1B. All the operation was performed in 
Discovery Studio 3.0. 


Molecular dynamics simulation 


All the molecular dynamics simulations were carried out with 
Amber 11, and the Amber 10 force field was used. Structures 
were solvated in a periodic truncated octahedral water box. 
The minimum distance between the complex and the box wall 
was 10 A. The TIP3P water model was applied and sodium 
ions were placed to neutralize the system. Bonds linking 
hydrogen atoms were restrained by the SHAKE algorithm. The 
particle mesh Ewald (PME) method was used to calculate 
the electrostatic interaction. Minimization of the systems con- 
sisted of two stages: at first, 5000 steps of steepest descent and 
conjugate gradient, respectively, were performed to avoid 
contact between water molecules by constraining the Ago/ 
siRNA complex. Then, the whole system was minimized 
under 5000 steps of steepest descent and conjugate gradient. 
After minimization of the solvated system, the general heating- 
equilibrium-production pipeline was performed on each struc- 
ture. The heating process was carried out for 20 ps from 0 K to 
300 K with 10 kcal mol’ A~* harmonic constraint under the 
NVT ensemble. The following equilibrium process lasted 
for 100 ps when the harmonic constraint was reduced to 
0 gradually under the NPT ensemble. Finally, 24 ns production 
was conducted under the NPT ensemble and the final 20 ns 
trajectory was used for analysis. The integration time step was 
1 fs and the coordinate was saved every 10 ps. 


Principal component analysis 


To observe the dominant motions of functional domains of 
Ago protein, principal component analysis (PCA)”° was carried 
out using AmberTools 1.5. First, the coordinate covariance 
matrix of Ago protein was built up for each complex. Then, the 
eigenvectors and eigenvalues of the covariance matrix were cal- 
culated by diagonalizing the matrix. The eigenvectors rep- 
resent the collective motions of the Ago domains, whereas 
the corresponding eigenvalues suggest the importance of the 
motions. In this study, the three most important motions were 
selected and projected onto the trajectory. The extreme projec- 
tion of each motion was used to indicate the directions of the 
functional domains. 


MM/PBSA binding free energy calculation 


The binding free energy between a ligand and a receptor is 
calculated by the following equation: 


AG» s = AGpy AGs. (AGs1 AGssr) 


where AG, and AG,, are the binding free energies in 
solution and vacuum, whereas AGs, AG, and AGs, are the 
solvation free energies for the complex, ligand and receptor. 
AG», consists of interaction energy between the ligand and 
the receptor, whereas sometimes an entropy contribution is 
also taken into consideration. Solvation free energy is com- 
posed of an electrostatic part and a hydrophobic part. In this 
study, the binding free energies of RISC/ssRNA for Ago/siRNA 
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complexes were calculated using the MM/PBSA method by the 
MMPBSA.py module in AmberTools 1.5. Because the differ- 
ences between complexes are only the isonucleotides, the 


entropy contribution is neglected. 
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